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Reconstruction of a surface topography , 

H "'■ 2003 



The invention relates to a method of reconstructing a surface of an object, 
where the object is represented by a 2-dimensional grid of measurements. For each grid point 
the measurements include corresponding information on a first slope of the surface in a first 
direction and a second slope of the surface in a different second direction. The invention 
further relates to software for executing the method and to a system for reconstructing the 
surface of an object. 



Surfaces of objects, such as silicon wafers, can be reconstructed from 
measured slopes. The slopes are typically measured by scanning the surface along lines and 
measuring the slope. The measurements result in a 2-dimensional grid, where for each grid 
point a slope in a 'horizontal' direction and a slope in a Vertical' direction is given. The 
measurements may be obtained using deflectometry, where light (e.g. from a laser) is 
projected onto the surface and an angle of reflection is measured, providing information on 
the slope. Fig.1 shows a typical equidistant measurement grid with measurements points 
along the rectangular coordinate lines. 

For many applications the actual surface (i.e. the topography of the surface) 
needs to be determined based on the measured slopes. Current reconstruction algorithms are 
based on carrying out a line integral along a path through the measurement grid. For 
examptej-by performing such a line integral along each 'horizontals-path (ie. each path 
parallel to one of the coordinate lines) the topography can be reconstructed. For each point 
aaa ^ pa*?* 6 integral uses the slope measured at the grid point in the direction of 
the path hi general the measured slopes contain measurement errors, Integration along the 
path accumulates all measurements errors in the slopes in the direction of the path. In this 
way, the topography of the surface determined at a grid point is therefore not determined 
uniquely but depends on the path chosen though the grid. For example, a line integral along a 
closed path like the path from (ij) to (i+lj), to (i+l, j+l) to (i, j+1) and back to (i j) will in 
general not give the starting value for the surface topography. 
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2 13.01.2003 
In general, the measurements will not result in the regular equidistant 
measurement grid of Fig. 1 with measurements points along the rectangular coordinate lines. 
In practice, measurement scan lines are not always straight but may be curved. If a 
measurement system is used that obtains slope information only one direction during a scan, 
scanning needs to be performed in both directions. In such a case, measurement points in one 
direction may not be perfectly aligned with the measurement points in the other direction. To 
compensate for irregularities in the measurements, usually a calibration procedure is 
performed to correct the position of the measurement points to obtain the perfect equidistant 
rectangular measurement grid. Such calibration is usually based on interpolation that in itself 
introduces further errors in the measurements and is time-consuming. 



It is an object of the invention to provide an improved method of 
reconstruction of a topography of a surface. It is a further object to provide software for 
performing such a method and a system for executing the method. 

To meet the object of the invention, the method of reconstructing a surface of 
an object, where the object is represented by a 2-dimensional grid of measurements, where 
for each grid point the measurements include corresponding information on a first slope of 
the surface in a first direction and a second slope of the surface in a different second 
direction, includes selecting a 2-dimensional part of the grid and fitting a corresponding part 
of the surface to the measurements of all grid points in the selected part, where the fitting for 
each grid point of the selected part is based on both the corresponding first and second slope 
information. By performing a fitting based on all measured slopes of the selected part, the 
effect of a localized measurement error is significantly reduced to mainly the area of the 
error. According to thememod, much more measurements are fully used for the - 
reconstruction. Instead of performing an integration operation along a one-dimensional path, 
now a 2-dimensional fitting operation is performed, involving much more measured data. 
Moreover, instead of using only one measured slope of each involved grid point, now both 
slopes are used. This enables significant reduction in propagation of measurement errors. 
Preferably, a user can indicate the area in which accurate reconstruction according to the 
invention is desired. Outside the selected area, a conventional reconstruction may be carried 
out. In a preferred embodiment, as described in the dependent claim 4, the fitting according 
to the invention is performed over substantially the entire surface represented by the grid 
measurements. 
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According to the measure of the dependent claim 2, the fitting is performed 
through a least-square minimization operation. This is an effective way of minimizing fitting 

errors. 

According to the measure of the dependent claim 3, the least square 
minimization operation is performed by solving an equation that describes a shape of a soap 
film loaded with a pressure field equal to a divergence of a slope vector including the first 
and second slope information. The inventor had the insight that in this way the surface fitting 
according to the invention can be expressed in a way similar to describing a shape of a soap 
film loaded with a pressure field. This enables using known methods for determining such a 
soap film shape to determine the topography of the surface. 

According to the measure of the dependent claim 5, for each point of the grid 
the first and second slope are measured using deflectometry. 

To meet the object of the invention, a computer program product operative to 
cause a processor to perform the steps of the method as claimed in claim 1. 

To meet the object of the invention, a system for reconstructing a surface of an 
object includes an input for receiving a 2-dimensional grid of measurements representing a 
surface of an object, where for each grid point the measurements include corresponding 
information on a first slope of the surface in a first direction and a second slope of the surface 
in a different second direction; a processor for, under control of a program, selecting a 2- 
dimensional part of the grid and fitting a corresponding part of the surface to the 
measurements of all grid points in the selected part, where the fitting for each grid point of 
the selected part is based on both the corresponding first and second slope information; and 
an output for providing a representation of at least the reconstructed surface part. 

According to the measure of the dependent claim 8, the system includes a 

measurement unit for measuring for each measurement point of a measurement grid the 

corresponding first and second slope information. 

According to the measure of the dep endent claim % the measuring is 
performed along non-straight lines; the measurement grid being directly used for the 
reconstruction. The method and system according to the invention perform a fitting that 
works as long as a connectivity between measurement points (analogous to the nodes 
belonging to a finite element in a FEM mesh) can be established. It is not required that the 
grid is neatly arranged, e.g. that the measurement points in the two directions form an 
equidistant x, y or r, <j> grid). No calibration of the measurement grid to a grid used for the 
reconstruction is required, avoiding the introduction of additional errors. 
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According to the measure of the dependent claim 10, the measurement unit 

includes a deflectometry measurement unit. 

These and other aspects of the invention are apparent from, and will be 

elucidated with reference to the embodiments described hereinafter. 



In the drawings: 

Fig. 1 shows a quadrilateral finite element mesh corresponding to the 
measurement grid; 

10 pig.2 shows ablockdiagram of asystemin which the invention can be 

employed; 

Fig.3 shows a finite element mesh consisting of triangular elements; 
Fig.4 illustrates a piecewise bilinear base function on element mesh shown on 
the 4 elements; 

15 Figs.5A and B shown an exemplary initial surface and reconstructed surface 

obtained from exact slopes with the method according to the invention; 

Fig.6 shows the discretization error, i.e. difference between reconstructed and 

exact prescribed surface; and 

Fig.7 shows the exact slopes, slopes containing noise used to reconstruct the 
20 surface, reconstructed surface and error in reconstructed surface. 



' Fig.2 shows a block diagram of a system in which the invention can be 
employed. In this exemplary system, the reconstruction of the surface topography is 

25 - performed by a processor 100 that is loaded with-a-suitable program. The program maybe 
stored in a permanent storage 150, such as a hard disk, and be executed from random access 
memory, such as DRAM (not shown). The processor operates on measurements received 
through an input 150. The input may be received from a measuring system 130. The 
measuring system 130 may be external to the system in which case the measurements may be 

30 received through a suitable communication system (e.g. wired or wireless LAN or wide area 
network, such as Internet). Preferably, the measuring system is part of the system for 
reconstruction of the surface. Advantageously, the measuring system is a deflectometry 
system, for example like the one described in US or US hereby included by reference. For the 
description it is assumed that measurements are provided for a 2-dimensional grid of 
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measurements points, where for each point a slope in the x-direction and the y-directioh is 
given. The method works equally well for other suitable coordinate systems, e.g. using an r, 
4> grid. It is not required that the grid is neatly arranged, e.g. that the measurement points iii 
the two directions form an equidistant x, y grid; it is sufficient that connectivity can be 
established between measurement points (analogous to the nodes belonging to a finite 
element in a FEM mesh). 

The input 120 is also able to receive input from a user. Shown are input 
devices, such as a mouse 140 and keyboard 150. Output to a user may be supplied via an 
output device 160, that may include a display. In particular the input may enable the user to 
indicate an area of the surface (and/or grid) in which accurate reconstruction of the surface is 
required. For other areas, conventional reconstruction may take place. Such conventional 
reconstruction may be used to reconstruct the edge of the selected area. In the remainder, the 
reconstruction method according to the invention will be described. 

Reconstruction method: ( 

Consider a domain D c « 2 on which the topography of an unknown surface 
is measured with a deflectometry measurement system. Using this measurement system, the 
slopes of the unknown surface z = z(x) are measured. If no measurement errors would 
occur, me unknown function mat describes toe surface satisfies 

#(*X XeD (1) 

Hence, the surface z(x) is determined apart from a constant function z = z 0 . A unique 
surface is obtained by prescribing the surface at some point X 0 eZ>,hence z(x 0 )=z 0 . 

In general the measured slopes contain measurement errors. Then, it is most 
likely that no function z(x) exists that satisfies(l), because in all cases where 

rot(ft)=V x ivV 0 no function exists that satisfies (1). However, it is possible to integrate 

(1) along a path r(x,x 0 ) from x 0 to X in the domain D to obtain values for the surface 
topography 

z(x)=z 0 + pv v -c& (2) 
rM 

Note, that if rot(ft)*0 these values are dependent on the path r(x,X 0 ) and are therefore not 
unique. This is seen by carrying out the line integration (2) along a closed path r(x 0 , X 0 ) 
enclosing a surface Q r and using Stokes theorem hence 
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J n dS=\\rot{$)-n Q da < 3 > 

where ft a is the unit normal on Cl T in tt 3 . Next, a method is proposed to find a unique 
function that approximately satisfies (1). Hence, the reconstructed surface obtained by the 
proposed method is independent of the fact that rot(fi) is zero or not. 

The unknown surface topography is approximated by a function z(x) that 
5 satisfies (1) in a least squares sense, hence it is the function with Z (* 0 )= z 0 that minimizes 
the functional 

j{z)=pz-ti)-$z-ri)da W 

D 

where the integration is carried out over the measurement domain D . In order to be a 
minimum for every disturbance Sz{x) of z{x) the minimizing function z{x) has to satisfy 

D 

This last formulation of the solution of the minimization problem (4) is the weak formulation 
10 of 

V.?z = V-#, XeD 
with boundary condition 

z(xo)=zo> n-Vz^n-R, xedD CT> 
where ft is the unit normal at the boundary 3D . 

The formulation (6) is also directly found by taking the divergence of the left 
and right hand side of (1). The equation (6) also describes the solution of the mechanical 
15 problem of determining the shape of a soap film loaded with a pressure distribution equal to 
V • n . Note, that in the weak formulation (5) the order of differentiation is one, while it is 
two in thelbriulation (6). Hence, solving (5) instead of (6) means that it is sufficient that the 
unknown.function *(*) Ues in the Sobolev space H x <p)- &D + *\ andis 

not necessarily in C x (d) (first order derivatives exist and are continuous). Furthermore, in 
20 the formulation (6) the first order derivatives of the measured slope # must exist. 

The advantage of solving (5) instead of carrying out 1he path integration (2) is 

that: 

- Both measured slope signals are used simultaneously to determine the 
topography of the measured surface. 
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- All measured slope data is used to reconstruct the topography of the 
measured surface. 

- The determined topography is independent of an integration path which has 
to be chosen using the simple line integration methods. 

- Analogous to the loaded soap film problem, the effect of errors 
(inaccuracies) in the measured slope data on the determined topography is mainly limited to 
the location of the inaccuracy. With the line integration methods a local error propagates 
along the remaining integration path and hence has a larger effect (larger error) on the 
reconstructed topography. 

Preferably, equation (5) is solved using a Finite Element Method (FEM). Of 
course other methods may be used to solve that equation, e.g. using a spectral method. 
Alternatively, equation (6) can be solved with e.g. a finite difference method. The 
formulation of the Finite Element Method to solve (5) will be carried out in Cartesian 
coordinates. For other coordinate systems the formulation is analogous to the presently . 
discussed one. 

Consider a Cartesian coordinate system (x,y,z). Suppose mat the 
measurement grid D u cDcdi 2 consists of n • m discrete points {x (l ,y JI ),i=l-n,J = l~. m , 
with x tJ < x MJ and y u < y MJ . With these discrete measurement points a quadrilateral finite 
elementmesh is established by connecting the points ( WJW ), (x MJ »,y^) 

and (*,. y+1 ,y y+w ) for * = l-«-l and y-l-w-l, as shown in Fig.l. Hence, the measurement 
domain is divided in (n-i).( m -l) quadrilateral elements D e using the nodes ( W J of the 
measurement grid. The choice of dividing the measurement domain D into quadrilateral 
elements is not a restriction. For instance, if each quadrilateral element is divided in two 
triangles, as-shown in Fig.3, a finite element mesh consisting of triangular elements can be 
created. In that case the number of elements increases with a factor 2. 

The unknown surface topography z = z(x,y) is written as aJinear combination 

of a finite set of approximating functions <p k (x,y) k = l—n-m, hence 

Substitution of (8) in (3) and niinimizing the functional j(z) with respect to the coefficients 
a k , or equivalently substitution of (8) in (4) with Sz = yields 
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(9) 

This singular set of linear equations is made regular by replacing the l* equation, where / is 
the index at which \ 9l (x 0 ) has a maximum, by the condition ^ 0 ) = Zo . The solution of this 
set of linear equations determines the approximation (8) of the measured surface topography. 
In the finite element approximation the subdivision of the measurement 
5 domain D in the finite elements D e is used to write the integrals in (9) as 

The functions 9t are chosen to be bi-linear per element and have the value 1 in one node and 
0 in all other nodes, as shown in Fig.4. Fig.4 illustrates a piecewise bilinear base function on 
elementmesh shown on the4 elements: [-l,o]x[-l,o], [o,l]x[-l,o], [o,l]x[o,l] and 

[-l,o]x[o,l]. 

10 For the evaluation of the integrals over the elements an iso-parametric 

transformation of an actual quadrilateral element to the standard quadrilateral element is 
carried out 

4 /-\ (ID 

where x = {x,yf , f = favf and the node number function is such that the actual 4 
nodes of an element are used. Then an integral of a function /(*) over an element D e is seen 
15 to be 



(12) 



where is the determinant of the Jacobian matrix F = J| . 
The bilinear approximating or basis functions <p k (4,v) seen to be 

(13) 

^,v>&+m-ri) (14) 

f,<fcifM*+f)G+*) (15) 
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d4 









(17) 




dtj 




(18) 




drj 




(19) 


^(1+^7); 


drj 


=i(w) 


(20) 



d4 

In the evaluation of the integrals in the right hand side of (9), integrals involving measured 
slope values occur. Using (10), (1 1) and (12) these integrals are transformed to integrals over 
the standard quadrilateral element. On the standard element the variation of the slope vector 
5 is approximated with a bilinear function 



(21) 



Substitution of (21) in the expression obtained by (10), (1 1) and (12) yields integrals of the 
form 

)WA^) d -^A4.n)mr, (22) 

-j -i og 



-1-1 or f 

These integrals can be evaluated analytically. However, in the present work a 4 point 
10 Gaussian quadrature rule is applied. The integration points and weights of this method i 

(67),=^V3,f^-^,=l (24) 

(£9)2 =4V3,f VJ); w 2 =l (25) 

=(h/3,*^)rs =1 (26) 

<g>v\=l&S,$j3}. w 4 =l (27) 
Carrying out all integrations, the element matrices and vectors are assembled 
in the large matrix and vector. Now the condition z(x 0 )=z 0 is incorporated into the 
equations, as discussed before and the resulting set of linear discrete equations 
4*«*=/, (28) 



are 
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has to be solved. Several techniques exist to solve (28). It is preferred to use a Gaussian 
elimination procedure (LU decomposition), which is a direct method. If the system of 
equations is very large, iterative procedures canbeused to solve (28), e.g. Gauss-Seidel, 

Multigrid, etc.. 

Examples: 

hi this section the accuracy and the efficiency of the proposed method is 
investigated using virtual measurement data; 

Suppose the measurement domain is D = [0,l]x [o,l]<= * 2 . The test surface is 

z(x,v)=sm 2 (2^)sin 2 (2^v) (29) 

with exact Slope vector 

(%\ _ ( 4*sin(2tf*)cos(2;rx)sin 2 (2* v)"| C 30 > 
[%) U^sin 2 (2^x)sin(2^r v)cos(2^v)J 
The initial surface and reconstructed surface obtained from the exact slopes (30) with the 
presently proposed method using an equidistant measurement grid of 21x21 points are shown 
in Figs.5A and B, respectively. The discretization error, Le. difference between reconstructed 
and exact prescribed surface, is shown in Fig.6. The reconstructed surface approximates the 
initial surface very good. Next, the accuracy of the obtained reconstructed surface is 
investigated when the measurement grid is varied, and/or noise is added to the slope vector. 

First, the discretization error to reconstruct the prescribed test surface with a 
finite discrete mesh is investigated. The discretization error to reconstruct the prescribed test 
surface with a finite discrete mesh are shown in Table 1, where f/j^ is the L 2 norm 

and is the infinity (or max) norm 

VL-i&kA (32) 

The discretization error reduces quadratically withthe mesh size (i.e. atypical dimension of 
and element, e.g. the radius of the largest circle that fits in an element). This behavior is 
expected using the bi-linear elements. 
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n 


m 


IK ~*mmL 


Ik -*«!L, 


10 


10 


76,808e-3 


37.788e-3 


20 


20 


24.838e-3 


9.2565e-3 — 


40 


40 


6.1787e-3 


2.3027e-3 


80 


80 


1.5428e-3 


0.57495e-3 



Table 1 : Discretization error to reconstruct the prescribed test surface. 

Secondly, using the same mesh (measurement grid) white noise (random 
values) with an amplitude equal to the amplitude of the slope values is added to the slope 
vector 

where a = ™*QrM- and r x (x) and r y (x) generate random values in the range 
In this case the difference between the reconstructed and the initial test surface is 

caused byadlscretizationerrorandbymenoisemmeslopedata. Fig.7 shows: Original . 
slopes: (A) in x-direction, (B) in y-direction; used slopes to reconstruct surface: (C) slope 
with noise in x-direction, (D) slope with noise in y-direction; and reconstructed surface (E) 
and error in reconstructed surface (F), for a grid consisting of 21x21 points. Even using the 
slopes with the added white noise having the same amplitude as that of the original slopes, 
the reconstructed surface approximates the original surface fairly good. As the problem is 
linear in the known slopes, the difference caused by the noise in the slope vector is linear 
with the amplitude of the noise. The difference between the reconstructed surface and the 
initial test surface is shown in Table 2. As follows from Table 1 and Table 2, the difference 
between the reconstructed and original surfaces is mainly caused by the noise in the slopes. If 
. the amplitude of the noise is halved, also the difference betweenjhe reconstructed and 
original surface halves. Note, that the behavior of the noise in the two considered cases of 
different amplitudes is different, because of the not repeatable random generated values. 
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10 



15 





noise amplitude = a 


noise amplitut 


ie = 0.5a 


n 


m 


IK -*«*L 




IK L 




10 


10 


0.69997 


0.24197 


0.28433 


0.08341 


20 


20 


0.41839 


0.12651 


0216 


0.061669 


40 


40 


0.31338 


0.075538 


0.13551 


0.042426 


80 


80 


0.18825 


0.050282 


0.084235 


0.018521 



Table 2: 



Difference between reconstructed and initial test surface, when noise is added to 



slope vectors. 

It should be noted that the above-mentioned embodiments illustrate rather than 
limit the invention, and that those skilled in the art will be able to design many alternative 
embodiments without departing from the scope of the appended claims. In the claims, any 
reference signs placed between parentheses shall not be construed as limiting the claim. The 
words "comprising" and "including" do not exclude the presence of other elements or steps 
than those listed in a claim. The invention can be implemented by means of hardware 
comprising several distinct elements, and by means of a suitably programmed computer. 
Where the system/device/apparatus claims enumerate several means, several of these means 
can be embodied by one and the same item of hardware. The computer program product may 
be stored/distributed on a suitable medium, such as optical storage, but may also be 
distributed in other forms, such as being distributed via the Internet or wireless 
telecommunication systems. 
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A method of reconstructing a surface of an object; the object being represented 
by a 2-dimensional grid of measurements, where for each grid point the measurements 
include corresponding information on a first slope of the surface in a first direction and a 
second slope of the surface in a different second direction; the method including selecting a 
2-dimensional part of the grid and fitting a corresponding part of the surface to the 
measurements of all grid points in the selected part, where the fitting for each grid point of 
the selected part is based on both the corresponding first and second slope information; 

2 - A method as claimed in claim 1, including performing the fitting through a 
least-square minimization operation. .-; 

3 - A method as claimed in claim 2, including performing the least square 
minimization operation by solving an equation mat describes a shape of a soap film loaded 
with a pressure field equal to a divergence of a slope vector including the first and second 
slope information. 

4 - A method as claimed in claim 1 , wherein the selected part of the grid is 
substantially the entire grid. 

5 - A method as claimed in claim 1 , including measuring for each point of the grid 
the first and second slope using deflectometry. ~ 

6 - A computer program product operative to cause a processor to perform the 
steps of the method as claimed in any one of the preceding claims. 



7 - A system for reconstructing a surface of an object including: 

- an input for receiving a 2-dimensional grid of measurements representing a 
surface of an object, where for each grid point the measurements include corresponding 
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information on a first slope of the surface in a first direction and a second slope of the surface 

in a different second direction; 

- a processor for, under control of a program, selecting a 2-dimensional part of 
the grid and fitting a corresponding part of the surface to the measurements of all grid points 

5 in the selected part, where the fitting for each grid point of the selected part is based on both 
the corresponding first and second slope information; and 

- an output for providing a representation of at least the reconstructed surface 

part. , 

10 8 A system as claimed in claim 7> wherein the system includes a measurement 

unit for measuring for each measurement point of a measurement grid the corresponding first 
and second slope information. 

9. A system as claimed in claim 8, wherein the measuring is performed along 
15 non-straight lines; the measurement grid being directly used for the reconstruction. 

10. A system as claimed in claim 8 or 9, wherein the system the measurement unit 
includes a deflectometry measurement unit. 
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ABSTRACT: 



A system for reconstructing a surface of an object includes an input 120 for 
receiving a 2 -dimensional grid of measurements representing a surface of an object. For each 
grid point the measurements include corresponding information on the slope of the surface in 
the x-direction and y-direction. A processor 100, operated under control of a program, selects 
5 a 2-dimensional part of the grid and fits a corresponding part of the surface to the 

measurements of all grid points in the selected part For each grid point of the selected part, 
the fitting is based on both the corresponding first and second slope information, An output 
160 of the system provides a representation of at least the reconstructed surface part to a user. 

10 Fig. 2 
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